A New mixed method with finite difference and finite element method for neutron diffusion calculation

内藤 俶孝; 鶴田 新一郎*; 林 正俊*

Journal of Nuclear Science and Technology, 18(8), p.571 - 580, 1981/00

https://doi.org/10.3327/jnst.18.571

一般に原子炉の幾何学形状はX-Y面は複雑であるがZ軸方向は単純である。このことに着目して炉内の3次元中性子束分布を得る新しい方法を開発した。この手法においては、X-Y面の計算には有限要素法が使用されZ軸方向には有限階差法が用いられている。3次元の中性子拡散方程式を解くのにこれ等2つの手法が交互に使用され、漏係数が矛盾なく求まるまで繰返される。この手法は計算機FACOM-M200用の計算コードFEDMとしてプログラム化されている。このプログラムを用いて3次元拡散計算が行なわれ在来の標準的な計算コードによる数値計算結果と比較された。その結果、両者は良好な一致を示した。また、本手法の有用性を示すために在来の階差法では計算困難な体系について計算も行なわれた。